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ABSTRACT 

The motion of a test particle in a stationary axisymmetric gravitational field is gen- 
erally nonintegrable unless, in addition to the energy and angular momentum about 
the symmetry axis, an extra nontrivial constant of motion exists. We use a direct 
approach to systematically search for a nontrivial constant of motion polynomial in 
the momenta. 

By solving a set of quadratic integrability conditions, we establish the existence 
and uniqueness of the family of stationary axisymmetric Newtonian potentials ad- 
mitting a nontrivial constant quadratic in the momenta. Although such constants 
do not arise from a group of diffeomorphisms, they are Noether-related to symme- 
tries of the action and associated with irreducible rank-2 Killing-Stackel tensors. The 
multipole moments of this class of potentials satisfy a "no-hair" recursion relation 
■^21+2 = a^M2i and the associated quadratic constant is the Newtonian analogue of 
the Carter constant in a Kerr-de Sitter spacetime. 

We further explore the possibility of invariants quartic in the momenta associated 
with rank-4 Killing-Stackel tensors and derive a new set of quartic integrability con- 
ditions. We show that a subset of the quartic integrability conditions are satisfied by 
potentials whose even multipole moments satisfy a generalized "no-hair" recursion re- 
lation M21+4 = (a^ + b'^)A{2i+2 —a^b'^M2i. However, the full set of quartic integrability 
conditions cannot be satisfied nontrivially by any Newtonian stationary axisymmetric 
vacuum potential. We thus establish the nonexistence of irreducible invariants quartic 
in the momenta for motion in such potentials. 

Key words: gravitation - black hole physics - celestial mechanics - chaos - stellar 
dynamics - galaxies: star clusters 



1 INTRODUCTION 

Since the work of Euler, Lagrange and Jacobi, the motion of 
a test particle in a Newtonian dipole fielcQ has been known 
to be completely integrable in terms of quadratures. In ad- 
dition to energy and angular momentum about the symme- 
try axis, there exists a third nontrivial constant of motion 
quadra tic in t he moment a. The constant was first discov- 
ered bv lEuleJ (|l760l . [1764I ) for two-dimensional (meridional) 
motion, and the problem is known as t he Euler problem 
l|Lukvanov. Emelianov. fc Shirrniiil 120051 ). Ibagrangel (117661 ) 
extended Euler's solution to three-dimensional motion and 
made a further generalization by allowing a Hookian (spring) 



* E-mail: charalampos.mark akis@uni-jena . de 
^ Following terminology in iMisner et al.l lll973l ). a Newtonian 
dipole will refer to a pair of fixed positive-mass centers each of 
which creates a Newtonian gravitational field 



center to be included between the two Newtonian (gravita- 
tional) cent ers. This generalization is known as the Lagrange 
problem (cf. Lukvanov et al.|[2005l ). The problem was further 
studied bv lJacobil l|2009l ) using separation of variables of the 
Hamilton-Jacobi equation in prolate spheroidal coordinates. 

The literature r egarding the problem o f two fixed cen- 
ters is survey ed by Lukvanov et al.l (|2005l ) and the orbits 
are studied in 16'Mathuna ( 20081 ). The quadratic^ constant 
has been studied by several authors using the Hamilton- 
Jacobi approach ([ St ackcl 1890; W hittakcr 1989; EddinrtOT 
I1915I : lKuzminlll956l : lLvnden-Beli|[r962l : Ide Zeeuwlil985al lbll3). 
If one considers the distance between the two centers to 
be imaginary, then one obtains a potential separable in 
oblate spheroidal coordinates. This is known in satellite 
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geodesy as the Vinti potential and has been used to ap- 
proximate the Kr a vitational field around the o blate earth 
(|Vintilll960l . 11961 Il969l . Il97ll : I Vinti et al.|[l99i l. A further 
generalization is possible by relaxing equatorial plane sym- 
metry. In the oblate case, this is accomplished by consider- 
ing two Newtonian fixed centers with complex-conjugated 
masses located at constant imaginary distance. The re- 
sulting Darboux-Gredeaks potential has been used to ap- 
pro ximate the gravitational field around o t her oblate plan- 
ets (Aksenov. Grebenikov. fc Deminlll963l . iLukvanov et al.l 
I2OO5I). 

Lvnden- Beil (|2003h has provided a simple and elegant 
derivation of the quadratic constant in the Euler problem, 
by noting that the kinetic part of the constant is the dot 
product of the angular momenta about the two fixed cen- 
ters of attraction. He also provided a generalization of the 
constant for a class of potentials that satisfy a certain in- 
tegrability condition. It will be demonstrated in section 
that this class is essentially that of the Lagrange problem, 
that is, it amounts to the addition of a Hooke term to the 
potential. 

[Israeli (|l970l ) and iKered ^96^ have demonstrated 
that the dipole field of the Euler problem can be re- 
garded as the Newtonian analogue of the Kerr solution 
in general relativity. Following Misner's suggestion to seek 
analogues of the quadratic constant in Newtonian dipole 
fields. Carter discovered his quadratic constant of mo- 
tion i n Ker r spacetime by separation of varia bles (|Carteil 
1 19681 . Il977l : iMisner. Thorne. fc Whe"eleil Il973h. The ana l- 
ogy has been furthe r eluc idate d by " iLvnden Belli (|2003D . 
iFlanagan fc Hindered (I2OO7I) a nd lWiUl (|2009l ) 



Carted (|l968l . l2009bl . l2010l ) also discovered a generaliza- 
tion of the Kerr solution with a non-zero cosmological con- 
stant, describing a rotating black hole in four dimensional 
de Sitter (or anti de Sitter) backgrounds. Carter's quadratic 
constant of motion exists for this class of spacetimes as well 
and has been used to s olve for the orbits in term s of hy- 
pergeometric functions (|Kraniotisll2004 l2005l . 120111 ). It will 
be shown in section 12.61 that the quadratic constant of the 
Lagrange two-center problem is the Newtonian analogue of 
the Carter constant in a Kerr-de Sitter spacetime. 

The existence of the Carter constant is very useful for 
analyzing the orbit of a small black hole around a mas- 
sive black hole. Detecting gravitational waves from these 
extreme mass ratio inspirals is a prime goal of the pro- 
posed New Gravitational Wave Observat ory, eLISA/NGO 
IJennrich et al.ll201ll : Umaro et al.ll2012al lbh. The infiuence 
of gravitational radiation reaction on the evolution of the 
Carter constant has been used to obtain the gravitational 
waveforms from orbits around Kerr black holes. As demon- 
strated bv iRvar] (|l995l ). by measuring the mass, spin and 
at least one more non-trivial moment of the gravitational 
field of a black hole candidate via gravitational wave ob- 
servations, one can test the validity of the no-hair theo- 
rem. This has raised interest in parametrizing spacetimes 
with mutipolc moments that deviate from those of the 
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times more general than Kerr-de Sitter that ad mit su ch 
a constant is an op en question ("Brink' '2008a"b', '2010aQ 
2OIII: iMirshekari fc Will 2010 : Kruglikov & Matvccv 20 111; 
Lukes- Gerakopoulosi 2012l V 



The aim of this work is to answer the analogo us question 
in Ne wtonian gravity. We use a direct approach (|Hietarintal 
Il987l ) to systematically search for a nontrivial constant poly- 
nomial in the momenta. This approach consists of directly 
solving the Killing equations and certain integrability con- 
ditions. In section[2]we show that the constant in the Euler 
and Lagrange two-center problems is the unique quadratic 
constant for motion around a stationary axisymmetric mas- 
sive object with equatorial reflection symmetry. In section 
I2.7l this constant is shown to be Noether-related. In section[3] 
we consider the next natural generalization, a constant quar- 
tic in the momenta. Using the direct approach, we show that 
no such constant exists for stationary axisymmetric vacuum 
potentials in Newtonian gravity. Conclusions are given in 
section |4| 



2 INVARIANTS QUADRATIC IN THE 
MOMENTA 

2.1 Killing equation and integrability conditions 

The motion of a test particle in a Newtonian gravitational 
field $ is independent of the particle mass, so for simplicity 
we can set the latter equal to unity. This motion is described 
by an action functional 



S[x,p] 



dt[pix'' — H{x,p)] 



(1) 



where position and momentum pi are treated as indepen- 
dent variables, 



H{x,p) = -g'' {x)p,pj + ^(x) 



(2) 



is the Hamiltonian, g^-' is the inverse of the Euclidian metric 
gij in E'' and summation over repeated indices is implied. 
Indices are raised and lowered with this metric throughout 
the paper. Unless otherwise noted, we will be using Carte- 
sian coordinates, so that gij = Sij where Sij is the Kronecker 
delta. The Cartesian components of the canonical momen- 
tum Pi = gijX-' will then coincide with those of the kinetic 
momentum x'. 

Spatial translations and rotations in are generated 
by the vector fields 



Xk — ek 



ry fc j 

±X,i — £ jiX Gk 



(3) 



Such deviations generally lead to loss of the Carter con- 
stant. Whether there exist stationary axisymmetric space- 



where Bk ~ {x,y,z} are Cartesian basis vectors and Sijk 
is the Levi-Civita tensor. We are interested in orbits of a 
test particle in the vicinity of a stationary and axisymmet- 
ric massive object, an object whose gravitational potential 
satisfies 

= (4) 

= if'^dk^ = {xdy - yd^)^ = (5) 

Here, £,p denotes the Lie-derivative along the Killing vector 
field 



<p = Rz = xy — y X 



(6) 
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that generates rotations about the 2:-axis (the axis of symme- 
try). By virtue of Noether's theorem, invariance with respect 
to time translations and azimuthal rotations implies respec- 
tively conservation of the Hamiltonian H and the component 

= if^'Pi = Rlpi = xpy - yp^ (7) 

of angular momentum. Rotations about the x and y axes, 
generated by the vector fields Rx and Ry are not considered 
symmetries of the problem, so the components = Pi 
and Ly = Ry Pi of angular momentum are not generally con- 
served. For example, the potential around a rotating massive 
object (such as a star or planet) in hydrostationary equilib- 
rium, may be expected to be axisymmetric, but not spher- 
ically symmetric, due to the rotationally induced deforma- 
tion. Without further symmetries, H and L~ are the only 
independent integrals of motion and, since the motion of 
test particles is three-dimensional, the problem is in general 
non-integrable. Nevertheless, one may seek special types of 
stationary axisymmetric potentials that admit a third non- 
trivial constant of motion and are therefore integrable. 

If one seeks a constant of motion linear in the momenta, 
then one is quickly led to rotational or translational symme- 
tries as the only choices, which have been already exhausted 
as explained above. The next natural step is to seek poten- 
tials admitting a nontrivial constant of motion quadratic in 
the momenta, 

I{x,p) = K''{x)p,Pj+K{x), (8) 

where the symmetric tensor K^^ and the scalar K are func- 
tions of position. The above quantity is conserved iff it com- 
mutes with the Hamiltonian, in the sense of a vanishing 
Poisson bracket: 

^^sr m = - = n (Q^ 

dt ^' '-dxf'dpk dpkdx'' ^' 

Substituting the Hamiltonian ((2]) and the ansatz ((8| into the 
above Poisson bracket yields 

{/, H} = d''K''p,p,pk + (d'K - 2K'''dk^)pi (10) 

where dk ~ tAt . In order that this Poisson bracket vanish for 
all orbits, the following nec essary and sufficient conditions 
l|Boccaletti fc PucaccdT2003l ) must be satisfied: 

^(fe^y) ^ (11) 

a'K = 2K'''dk^, (12) 

where parentheses denote symmetrization over the enclosed 
indices; that is, A^^^^ = \{A^^ + A^^). Eq. (p]) is a Killing 
equation in Euclidian space; a solution K^-' to the above 
equations wi ll be referred to as a Killing-Stackel tensor 
(|Carteilll977l ). From Eq. we obtain a necessary, but 

not sufficient, integrability condition 

d''\K'^'' dk^) = ^d^'ff^K = (13) 

where square brackets denote antisymmetrization over the 
enclosed indices, that is yl'*^' = ^{A'^ - A^'). 

Our assumptions of stationarity and axisymmetry al- 
ready guarantee the existence of two independent solutions 
to the above set of equations. First, the metric itself is al- 
ready a Killing-Stackel tensor, because equations pifl -HlS p 



are satisfied by K^^ = ip'^ = ij*^ , K = $, and the as- 
sociated conserved quantity is simply the Hamiltonian ((2}. 
Secondly, the above equations are also satisfied by the re- 
ducible Killing-Stackel tensor K^^ — f^'p-' , with K — 0, 
implying conservation of the quantity L?. This second case 
is reducible to the linear invariant ^ associated with the 
axial Killing vector ip'. We now explore the possibility of 
a third independent solution that can render the problem 
integrable. 

The conditio ns (|1HI - (|13(I suggest a systematic method 
(|Hietarintalll987l ) for obtaining integrable potentials and the 
associated invariants of motion: 

(i) The tensor if'-' may be computed by solving Eq. (|ll|l 
subject to the symmetries of the problem. 

(ii) With this tensor known, the integrability condition 
USjl provides a key restriction on the Newtonian potential 
One may solve this condition (e.g. via a multipole method) 
to obtain a family of integrable potentials. 

(iii) Given a family of potentials $ that satisfy this in- 
tegrability condition, one may obtain the scalar function K 
by integrating the components of Eq. (|12p . 

With this method as the basis of our analysis, we pro- 
ceed to carry out the prescribed steps in more detail. 

2.2 Rank-two Killing-Stackel tensors 

One may straightforwardly solve Eq. (|ll|l by noticing that 
the solution must be polynomial in the Cartesian coordi- 
nates. This is easily seen in one or two dimensions (c.f. Ap- 
pendix [XJ and c an be g eneralized to arbitrary dimensions 
and tensor rank (jHorwood 2008). In three dimensions, Eq. 
(Illfl constitutes an overdetermined system of ten equations 
for the six indepen dent components of the symmetric tensor 
K^-' . As shown by iHorwoodI (|2008l ). the most general solu- 
tion to this system is a sum of symmetrized products of the 
translational and rotational vectors (|3]): 

K = K''X,(g)X, = A''X^ (g) Xj+2B'' Xi®Rj+C' Ri<»Rj 

(14) 

The components K^-' of the tensor K are polynomial of sec- 
ond order in the Cartesian coordinates x^ and can be written 

as 

K^' = X' + 2B''^'e'\ix' + C^"e\^e\,,x''x' (15) 

where A""^ = C'^ = C**^^ are symmetric 3x3 constant 

matrices and B^-' is a nonsymmetric 3x3 constant matrix. 
While a superficial counting results in 21 independent coef- 
ficients, the actual number of independent coefficients is 20. 
This is because the main diagonal elements of B^-' appear in 
K'^ only in the combinations B^^-B^!', Bv^-B'"', B^^-B""" 
and the sum of these three terms is zero. 

2.3 Isometries 

One may considerably reduce the number of unknown coef- 
ficients by imposing known symmetries of the action func- 
tional S on the orbital invariant I. 

(i) Stationarity corresponds to invariance of S under the 
group action of (R, +) which represents time translations 
t ^ t + 6t. This symmetry has already been taken into 
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account, since all quantities have no explicit dependence on 
time t and an additive term dl/dt has been set to zero in 
Eq. ®. 

(ii) Axisymmetry corresponds to invariance of 5* under 
the group action of SO (2) which represents infinitesimal ro- 
tations {x, y,z} ^ {x + y Sip, y — X Sip, z} about the z-axis. 
This symmetry is generated by the Killing vector ((6| and 
may be imposed by requiring 



d 



d_ 

' dx 



d 



' dy " dx ' dpy dpx 

With I given by Eq. ([8]), the above condition can be shown 
to be equivalent to the requirement that the functions K'^ 
and K remain unchanged under rotations about the z-axis: 



d 



7 = 



(16) 



£^K = ip^dkK = 



(17) 
(18) 

Evidently, the symmetries (|16fl , (|17|l , (|18|) are inherited from 
the Hamiltonian, which is seen by the fact that these equa- 
tions are also satisfied by replacing I, K^-' , K with H, , "1? 
in the above three equations respectively. Substituting the 
general solution (|15p into Eq. (|17[) yields the axisymmetry 
constraints 



.4^ 



- = B""" + = 
: c''" = c - = 







(19) 

(20) 
(21) 



(iii) Further simplification is possible by assuming equa- 
torial plane reflection symmetry. This corresponds to invari- 
ance of S under the discrete group action of Z2 which repre- 
sents refiections {x,y, z,px,Py,Pz} {x,y,-z,px,Py,-pz} 
about the equatorial pla ne. (This assumptio n is not neces- 
sary for integrability, c.f. iLvnden-Belll (|2003l ). but we retain 
it for simplicity). Imposing this symmetry on the invariant 
(151) leads to the constraints 







= 


(22) 








(23) 




= s^" 


= 


(24) 






= 


(25) 



Note that the constraints (IT9l)-(l2TJ and (l22}-(l25} are 
independent, since axisymmetry and refiection symmetry 
are separate assumptions. Imposing these two sets of con- 
straints on the tensor p5|) yields 



K 




This is the most general solution to Eq. ([TT} consis- 
tent with the symmetry (R, +) x S0(2) x Z2. For clar- 
it y, we set k = C^ ^, X = ^ = C" - C^^ and a = 

■y/ {A''^ — A'^y)/Cyy. The parameter a is allowed to be real or 



imaginary and will be shown to depend on the gravitational 
source. (This reparametrization entails no loss of generality, 
as the sign of {A^^ — A^^)/C^^ is unrestricted. We have also 
excluded the possibility of vanishing with nonvanishing 
(^^^ - A^^). This case would lead to a Killing-Stackel ten- 
sor {A'"' - A^'")5l5i, reducible to a Killing vector SI that 
generates translations along the 2-axis. But this is not a 
symmetry of the problem by assumption, i.e. we have ex- 
cluded cylindrical symmetry.) The above solution may then 
be written as 



I + fj. (fi (fi-' 



(27) 



where gij — Sij is the Euclidian metric, tp is the axial Killing 
vector 



A'- 



r 2 ij 
(r g 

mn i j ( - 
■ 9 ^ km^ In \^ ^ 



2 ij i i\ 

a g — a a ) 



XX-') 

mn 1 j / k I k I 



- a a ) 



mn i i / k , k \ / I l\ 

■ g £ km,£ + a ){x - a) 



(28) 



is a new nontrivial Killing-Stackel tensor, r = y/x^Xk is a 
radial coordinate and a' = aS], are the components of the 
vector a — az. With K,X,jj, regarded as arbitrary coeffi- 
cients, the general solution (|27|l is a linear combination of 
three independent solutions: the known Killing-Stackel ten- 
sors 7'-' = S^-' and ip^ip'' , associated with stationarity and ax- 
isymmetry, and a third independent solution A^-' associated 
with a Noether symmetry to be discussed later. The above 
expression gives the most general rank-two Killing tensor in 
stationary-axisymmetric vacuum potentials. This completes 
step (i) of the prescribed procedure. 

2.4 Integrability condition: no-hair relation 

Although the solution (|28|l satisfies Eq. (Illf) and the symme- 
tries of the problem, it does not lead to a conserved quantity 
unless the condition (|12|) is satisfied. Our next step is thus 
to use the integrability condition (|13|l of Eq. (|12|) to obtain 
a family of potentials <I> for which A'^ leads to a conserved 
quantity. The Newtonian gravitational field around an ax- 
isymmetric object is completely characterized by a set of 
mass multipole moments {Ml}, by means of the expansion 



L=0 



(29) 



where Pl are the Legendre polynomials. This expansion is 
consistent with a stationary axisymmetric vacuum potential 
that vanishes at infinity. 

If one also requires equatorial plane reflection symme- 
try, then the potential "1? must be an odd function of z, so 
the odd moments Mi, Afs, . . . must vanish. Then, substitut- 
ing /f'-' = A^^ and $, as given by eqs. ((28l) and ((29)) . into 
the integrability condition (jlSp . we find by straightforward 
algebra that the latter is satisfied if and only if the even 
multipole moments satisfy the recursion relation: 

Ml+2 = o^Ml (L = 0, 2, . . .) (30) 

or, equivalently, 

Afi=a^Mo (L = 0,2,...). (31) 

where denotes the L-th power of a in the above two 
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equations. This relation is analogous to the "no-hair" re- 
lation for Kerr black holes (|lsraellll970l : lHansen|[l973 : IWilll 
I2OO9I ). All nonzero multipoles are determined by the mass 
Mo and the parameter a. If a is real (imaginary) then the 
quadrupole moment M2 = o^Mq is positive (negative) and 
the expansion (|29|l describes the field of a prolate (oblate) 
object. Summing the Legendre series (|29|l using Eq. (jSlf) 
yields (|Trahanaji20o3 ) 



Mo/2 



Mo/ 2 



+ + (z — a)^ + y'^ + {z + a)^ 



(32) 



In the prolate case, the above potential is that of Euler's 
three body problem: a test mass moving in the gravitational 
field of two point sources, each of mass Mo/2, fixed at posi- 
tions ±a z. 

In the oblate case, the replacement a — >■ ia (where i^ — 
— 1) allows one to write the above potential as 

Mo 



p + a'^z^/p'^ 

where p is an ellipsoidal coordinate defined by 



2,2 2 
X + y z 



= 1 



(33) 



(34) 



-I- p^ 

IVintil (|l960l . ll963l . [l969l , ll97ll ): IVinti et all (| 19981 ) has shown 
that the above potential is the most general solution to 
the three dimensional Laplace equation that seperates the 
Hamilton-Jacobi equation in oblate spheroidal coordinates. 
(A similar statement can be made for the potential 1)32^ in 
prolate spheroidal coordinates). Vinti's potential has been 
frequently used in satellite geodesy to approximate the grav- 
itational field around the oblate earth. 

We have so far shown that the solutions (|28|) and (|32p 
are the unique stationary, axisymmetric, equatorially sym- 
metric, vacuum solutions to eqs. (|lip and (|13fl . (Note that 
we have excluded the possibility of cylindrical symmetry as 
we are interested in the gravitational field of massive objects 
with compact support near the origin.) This completes step 
(ii) of our procedure. 

2.5 Existence and uniqueness of the quadratic 
invariant 

The last step towards obtaining an invariant is to obtain 
the scalar contribution K to the invariant ((8|. This is done 
by solving the condition p2p with A''^ = A^^ and $ given 
respectively by eqs. (|28|l and (|29p . The solution K — A is 
obtained by integrating the x component of Eq. (|12l) with 
respect to x, or the y component with respect to y. Up to 
some additive constant, we find 



A = 2 I A'^a,* = 2y dyA"'^^$ 

Moa(z + a) Moa{z - a) 



, ^ (35) 

x'^ + y"^ + [z + a)~ x"^ + y'^ + [z — aY 

Finally, substituting Eq. (|28|) into ((S]) and using the La- 
grange identity, the quadratic invariant is written 

/ = A'^p.Pj + A 

= t^IpI^ — \x ■ — (a^lpl" - |a ■ ) + A 

= \x X p\'^ — |a X + A 

= [{x + a) X p] ■ [{x - a) X p] + A (36) 



with A given by Eq. (|35|l . The above expression holds for 
the prolate case, a = \a\. In the oblate case, a = i|a|, 
the above expression still gives a true result. In the spheri- 
cally symmetric limit, a — >■ 0, the above invariant reduces to 
/ — >■ |x X p|^ which is a natural consequence of conserved net 
angular momentum. The first intergrals H,L^,I are inde- 
pendent and in involution (that is, {H, L^} ~ 0, {H, 7} = 
and {7,I/z} = 0). Thus, the system is Liouville-integrable. 

We have shown that the potential H33p is the unique sta- 
tionary, axisymmetric, equatorially symmetric, vacuum po- 
tential admitting a nontrivial constant of motion quadratic 
in the momenta, Eq. (|36p . The p otential (I33p can be con- 
sidered the Newtonian analogue (|Keresl Il967l : llsraell [l970l : 
lLvnden-Beiill2003l : IWilll 120091 ') of the Kerr solution in gen- 
eral relativity. The analogy is manifest when the Kerr met- 
ric i s written in Kerr-Schild coordinates (jKerr fc Schildl 
Il965h . The associated constant of motion is traditionally ob- 
tained by separating th e Hamilton-Jacobi equation in Bo yer- 
Lindquist coordinates (|Cartei|[l968l : iMisner et all 1 19731 ). It 
was in this way that Carter originally discovered his con- 
stant, following Misner's suggestion to seek analogies to the 
constant (|36p in Newtonian dipole fields. This analogy per- 
sists in the presence of a cosmological constant as demon- 
strated in the following section. 

Another interesting property is the following: One may 
interchange the roles of the metric gij and potential "I> with 
those of the Killing-Stackel tensor Aij and scalar field A 
respectively. Then the momentum map is generated by 7 = 
A^-' PiPj+A (which plays the role of the Hamiltonian) and the 
quadratic constant of motion is given by 77 = \g^''PiPj + 3>. 
This duality also exists in the relativistic case of the Kerr 
spacetime 



ty also exists ii 
|Carteij|2009"al ) 



2.6 Cosmological constant: analogy with a 
Kerr-de Sitter spacetime 

One may generalize the result of the previous section in 
various directions, by relaxing certain assumptions. The as- 
sumption that the motion is in vacuum may be relaxed by 
adding terms with a non- vanishing Laplacian to the multi- 
pole expansion (|29p of the gravitational potential. For ex- 
ample, one may attempt to add a spherically symmetric 
power-law potential of the form r'. Doing so, and repeat- 
ing the previous steps, leads to the same expression (|28p for 
the Killing tensor A'^-' as before. Direct substitution shows 
that the integrability condition (|13p is satisfied if and only 
if the moments {Ml} satisfy the no-hair relation (|30p and 
the only power law term that can be added is proportional 
to r^ . This leads to a potential 



Mo/2 



Mo/2 



x'^ + y^ + {z + a)2 x'^ + y'^ + {z — a)^ 
Aix^+y^ + z'') 



(37) 



In the prolate case, with a real, the above potential is that 
of two force centers, with Newton's inverse square law of 
gravitational attraction supplemented by a linear cosmolog- 
ical constant (or Hooke) term. (Note that the last term in 
the above equation may be regarded as the contribution of 
a spring connecting the test particle to the origin or as the 
contribution of two springs connecting the test particle to 
the two centers at ±az). The associated constant of motion 
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is given by Eq. (|36|) with 



A : 



Moa{z + a) 



Moa{z 



+ 4_ (2; + a)^ y/^;^ + + — a)^ 
+ 2Aa^(x-^ + y2) (38) 



In the oblate case, with a imaginary, the last term may 
again be interpreted as the contribution of a cosmological 
constant. With the substitution a — ia, the potential (|37p 
can be naturally regarded as the Newtonian analogue of a 
Kerr-de Sitter spacetime and the constant (|38|l as the ana - 
logue of the Carter constant in that spacetime ( CarteiifigGSl ). 
For this spacetime the metric has a tt component given by 
gtt = — (1 + 2$) in Kerr-Schild coordinates, with $ given by 
Eq. (I37II and a replaced by ia. 

iLvnden-Belll (|2003l ) provided a generalization of the 
quadratic constant p6p and the two-center potential (|32p . 
His result holds for potentials satisfying a certain integra- 
bility condition, which takes the form of a wave equation, 
{d^ /dr\ — /dr2){rir2^) = i n his two-center coordinates 
fi,2 = r ± a (|Lvnden-Bellll200^ . section 3). Assuming equa- 
torial symmetry and substituting the multipole expansion 
(|29p into this integrability condition, we recover the no-hair 
relation (|29p and the potential (|32p . If we attempt to add 
power law terms of the form ri , r2 to the potential $ we 
find that the only possibility is r\ -\- r2 ~ 2(r^ -I- a'^), giving 
rise to the potential (|37p up to a constant. We infer that 
Lynden-Bell's generalization accounts for the presence of a 
cosmological constant term, which gives rise to the New- 
tonian analogue of the Kerr-de Sitter spacetime discussed 
above. 

The direct approach employed here is more algorith- 
mic and computationally intensive co mpared to other ap- 
proaches (|Lvnden-Bellll2003l : I WilllbOOgI ) . But the advantages 
of the direct approach are that it establishes uniqueness and 
that it can be straightforwardly generalized to higher order 
or relativistic invariants. We have shown that the unique 
stationary, axisymmetric, equatorially symmetric potential 
admitting an invariant quadratic in the momenta is given by 
Eq. (I3 7p. The assumption of equatorial symmetry may be re- 
laxed (|Lvnden-Bellll2003l ') by dropping conditions (^^ - (^5)) . 
but we do not pursue this here. 



family of diffeomorphisms. Nevertheless, it is a generalized 
symmetry of the action and Noether-related to an invariant 
of the form gTJ. 

Conversely, if the quantity / is a constant of mo- 
tion, then the e- family of transformations generated 
by K''(x,x,t), obtained by solving the linear system 
(lloannou fc Apostolato3l2004l ') 



d^L ^ dl_ 

dx^dxi dxi ' 



(42) 



is a generalized symmetry of the action. 

The motion of a test particle in a Newtonian gravia- 
tional field can be obtained from the Lagrangian L{x,x) = 
^Qijx'x-' + If the motion admits a quadratic invariant 
I{x,x) = KijX^x^ -\- K, then the inverse Noether theorem 
(|42p implies that the e-family of transformations p9p gener- 
ated by = K^-'xj is a generalized symmetry of the action. 
We infer, for the problem of the previous section with A''^ 
given by Eq. (|27p , that the action ^ is invariant under three 
families of infinitesimal transformations: 

(i) The Killing-Stackel tensor ip^ip^ is Noether-related to 
the family of transformations 

x^ — X + eip if Xj = X + e Lz 

The tensor ip^ip^ is of course reducible to the axial Killing 
vector (p^, related to the family of diffeomorphisms xl = 
+ eif^. These represent azimuthal rotations and give rise 
to conservation of angular momentum ([7|. 

(ii) The metric g^-' is a Killing-Stackel tensor and is 
Noether-related to the family of transformations 

x^—x +eg Xj—x +ex 

or, equivalently, 

xi (t) ^ x\t) + € x\t) = x\t + e) 

The metric tensor g^-' is therefore Noether-related to invari- 
ance with respect to time translations t ^ t + e and gives 
rise to consevation of energy (or the Hamiltonian) given by 
Eq. ©. 

(iii) The irreducible Killing-Stackel tensor A^^ given by 
eq. (|28p is Noether-related to the family of transformations 



2.7 Generalized Noether symmetry 

The generalized Noether theo rem may be stated as follows 
l|loannou fc Apostolatosll2004h . Consider the e-family of in- 
finitesimal transformations 



xl = x^ + e K^{x, X, t) 



(39) 



which depend on position and velocity, for a small parameter 
e. If this family of transformations leaves the action S = 
J Ldt invariant, or, equivalently, changes the Lagrangian L 
by a total time derivative of some scalar K{x,t), 

dK 



L 



dt 



then the quantity 



dx^ 



(40) 



(41) 



is a constant of motion. Since the family (|39p of transforma- 
tions is velocity-dependent, it is not generally considered a 



or, equivalently, 

— X + e [{x X x) X X — (a X x) X a] 
= X + e [{x — a) X x] X {x + a) 

This a-posteriori knowledge of the symmetry transforma- 
tion allows a fast "derivation" of the quadratic invariant via 
the Noether procedure: varying the action of the two-center 
problem with respect to the above family of transformations 
yields no change, while the Lagrangian changes by —edA/dt 
with A given by eq. ((3HJ. Then, Eq. (|4ip leads immediately 
to the quadratic invariant (|36p . 

As mentioned above, these transformations are a symme- 
try of the action, giving rise to the constant of motion (|35p . 
but are not diffeomorphisms since they depend on position 
and velocity. Nevertheless, writing exj{t) = x-' {t + e) — x-' (t) 
they can be expressed as transformations in position and 
time: 

xi{t) = x\t) - A"xj(t) + A''xj{t + e) 
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where A^^ is a function of x'{t). 

The Noether theorem and its inverse, expressed by Eqs. (|39p - 
(|42p . also hold for a relativistic Lagrangian and can be used 
to show that the four constants of geodesic motion in a Kerr 
(or Kerr-de Sitter) spacetime are also Noether-related to 
symmetries of the action. Axisymmetry is related to con- 
servation of angular momentum about the symmetry axis 
as in case (i) above. Case (ii) discussed above has two ana- 
logues in general relativity: Stationarity (invariance of the 
Lagrangian under time translations along the integral curves 
of a timelike Killing vector) is related to conservation of en- 
ergy or Hamiltonian. Metric affinity (invariance of the action 
under proper time translations) is related to conservation of 
the magnitude of four-velocity (or the super-Hamiltonian) 
and is associated with the four-metric being a Killing ten- 
sor. Finally, a family of transformations analogous to those 
of case (iii) is relate d to the Carter constant of motion (c.f. 
|Padmanabhan|[201Cll . page 381). 



3 INVARIANTS QUARTIC IN THE 
MOMENTA 

3.1 Killing equation and integrability conditions 

The rank-two Killing-Stackel tensor (|28|) of the previous sec- 
tion reduces to a symmetrized combination of Killing vectors 
in the spherically symmetric limit a — 0. One might expect 
the next natural generalization to be a rank-four Killing- 
Stackel tensor, which would reduce to a combination of lower 
order tensors in some appropriate limit. We thus consider 
invariants quartic in the momenta 

Iix,p) = K'""{x)p,p,pkPi + K''{x)p,p, + K{x) (43) 

associated with a rank- four Killing-Stackel tensor A'*-**^'. The 
above quantity is conserved iff it commutes with the Hamil- 
tonian, in the sense of a vanishing Poisson bracket, Eq. @. 
Evaluating the bracket with I given by Eq. (|43|) and H given 
by Eq. ((U yields 

{I,H} = d'^K'^'^^piPjPkPiPm 

+ id^K^' - 4K'''''d,^)p,p,pk 

+ (a'=A--2/r'=9,$)pfc (44) 

Demanding strong integrability, that is, requiring that this 
Poisson bracket vanish for all orbits, we have 

Q(m,J^^jkl) ^ Q ^^^^ 

qC'K'') = 4K''''^di$ (46) 

d''K = 2K'''d^^ (47) 

As in the previous section, a solution to Eq. (|47|l exists only 
if the following integrability condition is satisfied: 

d^'iK'^'^dk^) = (48) 

which is identical t o condition (11311. T he above set of equa- 
tions is employed in lHietarintal |l983) to find quartic invari- 
ants for various two-dimensional systems. 

Eq. (|46p may be regarded as an inhomogeneous Killing 
equation and is also subject to an integrability condition. 
The general solution to this equation and its integrability 



condition are obtained in Appendix [X] We find that a solu- 
tion to Eq. (|46|) in exists only if the following integrability 
condition is satisfied: 

dyyyiK^'^'d,^) - 39„«. ( if^ "'9, $) (49) 

where i is summed over the y and z components and dij...k is 
an abbreviation for didj...dk- A generalization of the above 
condition to E" is straightforward and, as shown in Ap- 
pendix [X] given by 

^nml fijk Ojimifljk (^njl fijjik ~t~ dnjiflmk 

dkml fijn ~\~ (^kmifljn ~\~ dkjlfimn dkjiflmn — (50) 

where f^'' = AK^^'^^di^. Eq. (|49l) can be recovered from 
the above expression by setting n = m = i — z and i — 
j = k = y. Eq. (|49|l will suffice for our present purposes, 
as we shall restrict attention to two-dimensional motion in 
section |3^ and beyond. If, in addition to H and L^, a third 
constant of motion exists for all initial conditions in three 
dimensions, then this constant ought to also be conserved in 
the special case of orbits with = that lie on a meridional 
plane. That is, three dimensional motion is integrable only 
if meridional motion is integrable. We may thus set a; = 
and study motion restricted on the meridional (y-z) plane 
first. If such motion is integrable, generalization to three 
dimensions is straightforward by virtue of axisymmetry. 

Eqs. (|45|l - (|49|) suggest a systematic method for obtain- 
ing potentials admitting quartic invariants: 

(i) The tensor K^-'''' may be computed by solving Eq. (|45|l 
subject to the symmetries of the problem. 

(ii) With this tensor known, the integrability condition 
(|49|l provides a restriction on the Newtonian potential 
Solving this condition yields a family of possibly integrable 
potentials. 

(iii) Given a family of potentials $ which satisfy the con- 
dition (|49|l . one may obtain the tensor K''-' by solving the 
inhomogenous Killing equation (|46p . 

(iv) With K'^''^ and K'^ known, Eq. (gUl provides a fur- 
ther restriction on the potential ^. 

(v) If a family of potentials $ are found that satisfy 
the above integrability conditions, then one may obtain the 
scalar function K by integrating the components of Eq. (|47|l . 

We now proceed to carry out the above steps in more 
detail. 

3.2 Rank-four Killing-Stackel tensors 

In E'*, Eq. (|45|l constitutes an overdetermined system of 21 
equations for the 15 unknown independent components of 
the symmetric tensor K^^'"'- . Similarly, in E^, Eq. (|45)) is a 
system of 6 equations for 5 unknowns. The system may be 
solved by noticing that each component of K'-'''' must be 
a fourth order p olynomial in the Cartesian coordinates a;'. 
iHorwoodI (|2008l ) considered the equation a^^if = 
for a tensor ' Qf arbitrary valence A'^ and showed that 
the general solution is given by 

• ■ • (g) X,^ R,^^^ (gi ■ • ■ (gi R^j^ (51) 
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where the objects C^^ labelled by L = 0,...,N, are 
constant and subject to the symmetries C\ 

The above expression, for = 4 and 



(n---it){ii+i'"ijv) 



c 

with Xi, Rj given by Eq. ((3|, provides the general solu- 
tion to Eq. (|45|) . Imposing known isometries on the quartic 
invariant (|43p constrains the coefficients C^i "'™ 

Stationarity has already been imposed (as the system 
is autonomous). For three dimensional motion, axisymmetry 
may be imposed via the requirement that the above Killing- 
Stackel tensor is Lie-derived by the axial vector (p = R^, 
that is 



to 7 



- K 



mijl 



mip - K dm<P = 



(52) 



(For purely meridional motion, enforcing the above condi- 
tion does not change the relevant coefficients of K^-'''' pro- 
jected to the meridional plane.) In addition, we also re- 
quire li2 reflection symmetry about the equatorial plane, 
that is, the replacement {z,pz} — >■ {—z,—pz} leaves the 
quantity K^-'^^PiPjPkPi unchanged. Since the latter quan- 
tity is a fourth order polynomial in the (Cartesian) position 
and momentum variables, imposing the above isometries is 
straightforward via algebraic manipulation software such as 
Mathematica. Imposing the above symmetries signiflcantly 
reduces the number of independent polynomial coefficients. 
Then, after dropping trivial terms such as g^'-'p'^'', g^^^ ip^Lp^\ 
ip^ ifi^ tfi^ ip'' which correspond to reducible Killing-Stackel ten- 
sors, we find that the most general nontrivial solution to (|45|) 
subject to the symmetry (R, -f ) x SO(2) x Z2 can be written 
in the remarkably simple form 



K 



ijkl 



(53) 

where A'-* is given by Eq. (|28p and B*-* is given by the same 
equation with the replacement a — b: 



A'- 



mn i i / k , k \ / I l\ / 1- a\ 

J e km£ ini^ +a){x - a) (54) 



j-jii mn i i / k , j k \ / I j l\ 

B ' = g e km£ +b )[x - 6 ) 



(55) 



where a — a z,b — b z. Note that the above expressions are 
valid in two and three dimensions. This completes step (i) 
of our prescribed procedure. 

3.3 First integrability condition: no-hair relation 

With the tensor K''-'''^ known, the next step consists of find- 
ing the class of gravitational potentials $ that satisfy the 
integrability condition H49p . A stationary axisymmetric vac- 
uum potential that vanishes at infinity is characterized by 
the multipole expansion (|29|) . Assuming equatorial plane 
reflection symmetry, the odd moments Ml {L — 1,3,...) 
vanish. Restricting attention to purely meridional motion, 
with K'^'''- and $ given by eqs. (f55)) and we find with 
straightforward algebra that the integrability condition (|49p 
is satisfied if and only if the even multipole moments satisfy 
the two-step recursion relation 

Ml+4 = (a" + 6^)A/l+2 - a'b^Mi. (L = 0,2,...) (56) 

This equation generalizes the "no-hair" relation H30p . One 
may obtain all higher moments recursively from the first two 
nonvanishing moments. This yields 

(a^fe^ - a^b^)Mo + (a^ - b^)M2 



Ml = 



fo2 



(57) 



for L even (and Ml = for L odd). In the limit b — >■ 0, 
one recovers the no-hair relations (|30p and (|3ip (except for 
L = 0). The limit 6 — >■ a is slightly more subtle. Although 
the 6 — >■ a limit of eq. (|56p is satisfied by pO|) the converse 
is not necessarily true, because Eq. (|56p does not fix the 
relation between the first two moments. Applying Eq. 1)56^ 
recursively after taking the limit b a (or taking the same 
limit of Eq. (|57p and applying the 1' Hospital rule) yields 



lim Ml = a 



Mo 



^ia-'M2 



Mo) 



(L = 0,2,...) (58) 



If the mass and quadrupole moment are related by M2 = 
a^Mo then one recovers Eq. pip , but this need not be the 
case and M2 is generally considered independent of Mq. If 
a b, then the gravitational field depends on four indepen- 
dent parameters {Mo, M2,a,b) and the reparametrization 
rua = (M2 - bH4o)/{a^ - &^), mt = (M2 - aH'Io)/{b' - a') 
allows us to write Eq. (|57p in the suggestive form 

Ml = a^m, -f 6^7716 (L = 0,2,...) (59) 

Then, summing the Legendre series (|29p yields the potential 

ma/2 ma/2 



^x^ +y'2 + {z~ a)2 ^2,-2 +y^ + {z + a)2 

mb/2 mb/2 
^x^ +y'^ + {z- 6)2 ^x^ + {z + b)^ 



(60) 



created by four fixed point sources with mass ma/2 at po- 
sitions ±az and mass mb/2 at positions ±62. Our solu- 
tion procedure guarantees that the above expression gives 
the unique vacuum potential with the symmetry (R, -I-) x 
S0(2) X Z2 compatible with the integrability condition (|49p . 
This completes step (ii) of the prescribed procedure. 

3.4 Second integrability condition: nonexistence 
of quartic invariants 

With $ given by Eq. (|60}, one may use eqs. ()A6P - (|A9P to 
solve the inhomogeneous Killing equation (|46p . We find 

K'' =2BA" + 2AB'' 

+ {b^ - a^){A - B){z5i5i - ySi'Si^) + (61) 

where i/ is a constant, A is given by Eq. (|35p with the re- 
placement Mo — >■ 2ma, A is given by 

ma{z + a) ma{z ~ a) 



A = 



^y^ + {z + af ^y^ + iz^af 



(62) 



and B, B are given by Eqs. (|35p . (|62p with the replacements 
a — > fo and ma — >■ mb- The additive contribution C*-* to 
Eq. (|6ip is a solution to the homogenous Killing equation 
(jlip . which is polynomial of second order in the cartesian 
coordinates. Since the homogeneous solution must obey the 
symmetries of the problem, it must have the form of Eq. 
with a replaced by some other parameter c, that is 



C-'=g e km£ ini^ +c){x-c) 



(63) 



where c — cz. This completes step (iii) of our prescribed 
procedure. 

The next step is to solve the second integrability condi- 
tion (|48p . However, substituting Eqs. (|60p and (|6ip into the 
condition (|48p, we find that the latter cannot be satisfied 
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except in the limit b ^ a, whence the quartic invariant is 
reducible to the quadratic invariant of the previous section. 
We infer that there exists no stationary, axisymmetric, equa- 
torially symmetric vacuum Newtonian potential that admits 
an independent nontrivial invariant quartic in the momenta. 
That is, the only quartic invariants are trivial products of 
lower-order invariants. 

The proof of this nonexistence result uses the fact that 
strong integrability requires existence of a constant of mo- 
tion for all values of energy E and angular momentum L,, 
including the case of purely meridional orbits with = 0. 
Failure to satisfy Eq. H49p means that there exists no in- 
dependent constant for purely meridional motion; therefore 
there exists no strong integral for three-dimensional motion. 
One could alternatively use the integrability condition (|50p 
to derive the same nonexistence result directly in E''. 

The nonexistence of a quartic invariant does not pre- 
clude the possibility that the potential (|60|) admits some 
other constant of motion with different dependence in the 
momenta. A superficial study of Poincare maps of three di- 
mensional orbits in the four-center potential may show that 
most orbits appear to be regular. This could lead to the im- 
pression that the system is integrable. However, a thorough 
scan of initial conditions reveals the existence of Birkhoff 
chains and in some cases ergodic motion surrounds the main 
island of stability on the Poincare maps, confirming the non- 
integrability of the system described by the potential (|60p . 
In a relativistic context, a similar behaviour has been ob- 
served for orbits in certain stat i onary axisymmetric space- 
times (lApostolat os et al" 20091 ; iLukes-Gerakopoulos et al] 



l2010l : [Cantopoulos et al.ll201ll ') 



4 SUMMARY AND CONCLUSIONS 

The aim of this work was to study polynomial constants of 
motion in stationary axisymmetric gravitational fields in a 
Newtonian context and to explore the analogies with similar 
problems in relativistic gravity. The results in this study that 
are (to the author's knowledge) new include: 

(i) the uniqueness of the constant (|35l) - (|36| ) of the Euler 
problem: Eq. p2p gives the only stationary, axisymmetric, 
equatorially symmetric Newtonian vacuum potential admit- 
ting a Killing-Stackel tensor of rank two; 

(ii) the uniqueness of the constant ((36])- ((38]) of the La- 
grange problem: Eq. (|37p gives the only stationary, axisym- 
metric, equatorially symmetric Newtonian potential admit- 
ting a Killing-Stackel tensor of rank two; 

(iii) the relatio n of th e quadratic constant in the Lagrange 
problem to that of lLvndcn-Bc U ( 2003), and the analogy with 
the Carter constant in a Kerr-de Sitter spacetime; 

(iv) the integrability conditions (|49p . (|50p which do not 
appear to have been implement ed previously i n direct 
searches for quartic invariants (c.f. lHietarinta|[l987l ). 

(v) the non-existence of stationary, axisymmetric, equa- 
torially symmetric, vacuum Newtonian potentials admitting 
a Killing-Stackel tensor of rank four. 

Note that the assumptions of vacuum and equatorial sym- 
metry may be relaxed and th e two-center problem remains 
integrable (|Lvnden-BelJl2003l ). Note also that the integrabil- 
ity condition (|49p and its generalization to arbitrary dimen- 



sion, Eq. (|50l) . can be quite useful for other direct searches 
of quartic invariants. 

The methods employed here may provide intuition 
for treating similar problems in relativistic gravity or 
post-Newtonian approximations to it. This could lead to 
further insights on how to parametrize the departure of 
a spacetime geometry from that of a Kerr spacetime and 
the relation of this departure to integrabilty or ergodocity. 
Extrapolating the present results to the relativistic regime 
would lead to the conjecture that rank-4 Killing tensors do 
not exist in stationary axisymmetric vacuum space times. 
This conjecture would contradict claims made by iBrinkI 
(|201ll ) on existence of rank-4 Killing tensors in certain 
stationary axisymmetric spacetimes, but would be in agree- 
mcnt with the nonexistence results of Kruglikov fc MatveevI 



(2011) and Lukcs-Gcr akopoulosI (|2012l ') in the case of the 
Zipoy-Voorhees spacetime. Nevertheless, extending the 
Newtonian analysis to generic stationary axisymmetric 
vacuum spacetimes is a nontrivial step. 
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APPENDIX A: INTEGRABILITY CONDITIONS 
AND INHOMOGENEOUS KILLING 
EQUATIONS 

We consider the inhomogeneous Killing equation 

d(kK,j) = (Al) 

where P^'' = AK'^^'^' di$. In E^, the cartesian components of 
the above equation read 

9yKyy — fyyy (-^-2) 
dzKyy -\- IdyKyz ~ '^fyyz (^'^) 
dyKzz + 2dzKzy = ifyzz (A4) 

dzKzz = (A5) 

The above system is overdetermined and may be solved as 
follows: Integrating Eq. HA2|) with respect to y and Eq. HA5|I 
with respect to z yields 



Kyy = Z{z) + / dyf, 



Kzz = Y{y) + / dzf: 



(A6) 
(A7) 



where Z{z) and Y{y) are scalar functions of their arguments. 
Then, integrating Eq. (IA3P with respect to y and Eq. HA4|I 
with respect to z yields 



K. 



yz - C(2) + 2 / dyiZfyyz - dzKyy, 



(A8) 



Kz 



i>{y) + dziSfyzz - dyKzz) (A9) 

where C,{z) and a-re scalar functions of their arguments. 
Eqs. (|A6|) - (|A9|) provide the solution to the inhomogeneous 
Killing equation ([Aip . 

Because Kij is a symmetric tensor, the above two ex- 
pressions must be equal. Acting with dyyzz ~ Q^iQz'i on Eqs. 



[8|l . (|A9|) . demanding that the two expressions be equal 
and using Eqs. (|A6|I . (|A7|) . yields the integrability condition 

dzzzfyyy — Sdzzyfyyz = dyyyfzzz — Sdyyzfyzz (AlO) 
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which must necessarily be satisfied by fijk in order for 
Eq. (|A1[) to have a solution. Note that the unknown func- 
tions Z, Y, -0 do not appear in the above condition. In 
the homogeneous case, these functions can be easily shown 
to be quadratic polynomials in their arguments by setting 
fijk ~ 0, demanding that expressions (|A8I) and (|A9I) be 
equal and separating variables. 

The integrability condition ()A10[I may be generalized to 
E" as follows. Eq. (|AH) can be expanded out as 

diKjk + djK.u + dkKi, = 3/,,fc (All) 

The first term vanishes if we apply di and antisymmetrize 
over / and i. The second term vanishes if we subsequently 
apply dm and antisymmetrize over m and j. Finally, the 
third term vanishes if we apply dn and antisymmetrize over 
n and k. This yields the integrability condition 

dnml fijk f^nmi fljk (^njlfimk Onj i fhnk 

dkml fijn ~t~ dkrni fljn ^^ Okjl fimn dkjiflmn — (A12) 

which generalizes Eq. (|A10|) to arbitrary dimension. 

This paper has been typeset from a T^j]X/ ffl^jX file prepared 
by the author. 
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